Astronomy &: Astrophysics manuscript no. 14605 © ESO 2010 

May 28, 2010 



O 

(N 



Letter to the Editor 

In-flight calibration of the Herschel-SP\RE instrument* 

B. M. Swinyard 1 , P. Ade 2 , J-P. Baluteau 3 , H. Aussel 16 , M. J. Barlow 4 , G. J. Bendo 5 , D. Benielli 3 , J. Bock 6 , 

D. Brisbin 7 , A. Conley 8 , L. Conversi 9 , A. Dowell 1 , D. Dowell 6 , M. Ferlet 1 , T. Fulton 10 , J. Glenn 11 , A. 

Glauser 19 ' 20 , D. Griffin 1 , M. Griffin 2 , S. Guest 1 , P. Imhof 10 , K. Isaak 2 , S. Jones 14 , K. King 1 , S. Leeks 1 , L. 

Levenson 12 , T. L. Lim 1 , N. Lu 13 , G. Makiwa 14 , D. Naylor 14 , H. Nguyen 6 , S. Oliver 15 , P. Panuzzo 16 , A. 

Papageorgiou 2 , C. Pearson 1 ' 14 , M. Pohlen 2 , E. Polehampton 1 ' 14 , D. Pouliquen 3 , D. Rigopoulou 1,21 , S. 

Ronayette 1,16 , H. Roussel 17 , A. Rykala 2 , G. Savini 4 , B. Schulz 13 , A. Schwartz 13 , D. Shupe 13 , B. Sibthorpe 19 , 

S. Sidher 1 , A. J. Smith 15 , L. Spencer 2 , M. Trichas 5 , H. Triou 16 , I. Valtchanov 9 , R. Wesson 4 , A. Woodcraft 18 , 

C. K. Xu 13 , M. Zemcov 11 , and L. Zhang 13 

(Affiliations can be found after the references) 



J>^' Received May 28, 2010; accepted 

& . 

ABSTRACT 

*vj ' SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrome- 
ter. It contains a three-band imaging photometer operating at 250, 350 and 500 /im, and an imaging Fourier transform spectrometer 
i— i, (FTS) covering 194-671 /jm (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the 
^ ' SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission 
^H ' from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that 
™"3 \ is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission 
^S^ . is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the 
<**> i' photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The 
spectrometer wavelength accuracy is determined to be better than 1/10*'' of the line FWHM. The astrometric accuracy in SPIRE 
maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is 
currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing 
ffi ■ methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the 
instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models 
of astronomical standards. 
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c"^ ■ 1. Introduction Fourier transform spectrometer (FTS) covering 194-313/im 

IT) • , . „. , ,., . , . . . (SSW) and 303-671^m (SLW) with arrays of 37 and 19 

. ■ The m-flight calibration of any astronomical instrument re- detectors SPIRE contains tw0 internal ca ii bra tion sources 

^ lies on a combination of previous measurements using in- which were desi d to provide a rapidly modu i ated sig _ 

g . strumentation with known or traceable calibration accu- nal for rektive in calibration (the PCAL source) and a 

^ ; racy, accurate models of standard astronomical sources and gtable thermal SQurce for calibration of the spec trometer 

. calibration sources internal to the instrument or facility. In and tQ balance ^ pQwer from the telescope (the SCAL 

_> . this paper we discuss the methods ; and source models em- source) pcAL ig 1&ced ftt an { rf the tel e sec _ 

^ ; ployed to convert SPIRE data to physical units and discuss Qnd ^ ^ common Ught th of the tw0 instruments 

^ . the accuracy of the calibration and any caveats that must and therefore can be used to stimulate a ll detectors in the 

CG ■ be placed on the data. photometer and spectrometer. The SCAL source is placed 

The SPIRE instrument and its overall calibration ^ ^ second input pQrt of the spectrometer (Swinyard et 

scheme are described m detail in Griffin et al. tfMn). For &1 ^^ R ig therefore constant l y vie wed during observa- 

the purposes of discussion of the calibration and perfor- tiong and provides a re f er ence against which the spectrum 

mance estimation we treat SPIRE as two separate mstru- from ^ telescope port is mea sured. The intention before 

ments that are spatially separated at the focal plane of flight wag that the gcAL SQurce WQuld be heated tQ & tem _ 

the Herschel telescope: a three band imaging photometric perature sufficient to entirel « nuU » the spec trum from the 

camera with bands centred nominal ly on 250 M m (PSW), teles We discuss m Sect. |6] why this has proved to be 

350 T.iJ^'^ Mm n ( 1 e ^r^" unnecessary for standard observing conditions. The inter- 
rays of 139, 88 and 43 respectively feedhorn coupled NTD- na] calibrators are only viewed throu h part of the optics 

bolometers (Turner et al [MTJ , and a two band imaging chain and the Qverall astronomical ca i ibr ation and perfor- 

* TT , , . nn, , , .,, . • mance of the instrument and telescope together can only 

Herschel is an ESA space observatory with science mstru- . . 

ments provided by European-led Principal Investigator consor- be established by observation of astronomical sources as 

tia and with participation from NASA. w e describe in this paper. 
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The structure of the paper is as follows: in Sect. [5] we 
describe the initial in-flight performance evaluation of the 
instrument before going on to cover the basic conversion 
of the data from digital encoder values to physical units 
for the detectors, the beam steering mirror (BSM) and the 
spectrometer mechanism (SMEC) in Sect. [3] In Sect. 2] we 
describe the characterisation of the time response of the 
photometer bolometers using ionising radiation hits and in 
Sects. [5] and [B] we discuss how the photometric calibration 
was established for the photometer and spectrometer. In 
Sect.[7]we present the results of comparing SPIRE observa- 
tions to solar system objects with known or well constrained 
model flux densities. The bcamsize and astrometric accu- 
racy are discussed in Sect. [H] and the wavelength accuracy 
in Sect. [51 We summarise the present calibration status of 
SPIRE and draw some conclusions for future calibration 
activities in Sect. [TO] 

2. Initial performance evaluation 

The initial performance verification and comparison be- 
tween in-flight and ground performance was carried out us- 
ing the PCAL source whilst the Herschel cryostat lid was 
closed. This allowed a direct measurement of the detec- 
tor response under conditions similar to those seen dur- 
ing ground testing. The response measured using PCAL 
before opening of the cryostat lid showed a close corre- 
lation between ground and flight in all detector channels 
in both the photometer and spectrometer. The overall in- 
flight measurement showed a slightly higher response in all 
arrays as the bolometers are approximately 10 mK colder 
in flight than during spacecraft level ground tests. The time 
response of the detector channels is a combination of the 
thermal response of the bolometers and the low pass fil- 
ters in the readout electronics. In the photometer chan- 
nels the time response is dominated by the low pass filters 
which have a nominal 5 Hz 3 dB cut off. Here the PCAL 
flashes, which have switch-on time constant of 90 msec and 
a switch-off time constant of 50 msec, can be used to make 
a first order assessment of the response of the system and 
look for anomalous behaviour. This test is not possible for 
the faster (25 Hz) spectrometer filters and bolometers al- 
though analysis of the interferograms does provide a first 
order check. This is achieved by measuring the difference 
in recorded position of the peak in signal at zero path dif- 
ference when the spectrometer mechanism is scanned in 
its "forward" and "reverse" directions. In both the pho- 
tometer and spectrometer all channels that were working 
nominally before flight showed the same or similar time 
response. A more detailed investigation of the photometer 
time response was carried out using an analysis of ionising 
radiation glitches and is discussed in Sect. |H 

Once the lid was opened we were able to observe both 
astronomical sources with known or well constrained model 
flux densities and the emission from the Herschel telescope. 
The telescope primary and secondary mirror temperatures 
are measured using several thermistors with a notional 0.1% 
accuracy. This, in addition to measurements of the mirror 
emissivity (Fischer et al. 120040 . allows us to model the di- 
rect emission from the telescope which fills the field of view 
of all detectors. By varying the electrical power applied 
to the bolometers through changing the bias voltage - i.e. 
taking a so called "loadcurve" - it is possible to make a di- 
rect estimate of the optical power absorbed by the bolome- 



ter. Comparing this measurement to the expected absorbed 
power using the modelled and ground based measurement 
of the instrument throughput and transmission showed that 
a) there is little or no straylight in the SPIRE beam due to 
sources on the satellite and b) that the instrument through- 
put and transmission are as good or better than expected 
before flight in both the photometer and spectrometer. 

To assess the optical performance of the system and to 
calibrate the response of the photometer to point sources 
two separate types of scanned observations were carried 
out. The first takes a scan map observation using the stan- 
dard Astronomical Observation Template (AOT) of a point 
source of known or modelled flux density. In reconstructing 
the map from these data the time series from many differ- 
ent detectors in the array are combined to make a single 
image of the beam. The second approach was to scan a 
bright point source along the length of the array in a raster 
pattern with fine steps between the scan legs allowing the 
beam map, relative position and response of each individual 
bolometer to be determined. The results from both the pho- 
tometer and spectrometer arrays confirmed that the beam 
sizes are as expected and that no significant change to the 
optical alignment of the instrument had occurred during 
launch. These measurements are discussed in further detail 
in Sect. [8] 

3. Data processing and calibration 

The data processing pipeline for the SPIRE photometer 
and spectrometer data are discussed in detail in Griffin 
et al. (|20"0g|i and Fulton et al. (|20"0g|t . In this section we 
highlight those elements of the processing that require the 
application of calibration data and discuss the origin and 
accuracy of the applied conversion factors. 

The basic conversion from digitised output from the 
SPIRE electronics to voltage (in the case of the detectors) 
or position (in the case of the BSM and SMEC) relies on 
the accurate measurement of the gains and offsets of the 
amplification chains and encoders in the SPIRE warm elec- 
tronics. These were extensively tested during the ground 
calibration and again in flight. Fixed resistors mounted on 
the bolometer arrays were used to check the detector con- 
ditioning electronics gains and no change was found com- 
pared to the ground measurement. The values of the gains 
applied in pipeline processing can be assumed to much more 
accurate than the overall calibration based on astronomical 
flux densities and need not be considered in estimating the 
calibration accuracy. 

The SMEC position encoder accuracy was tested by 
comparing the measured and expected frequencies of known 
lines (see Sect. 0. The position encoder of the BSM has 
been recalibrated in flight by using it to move a bright point 
source across the field of view and using the known posi- 
tions of the detectors (see Sect . |HJ) to determine its position. 
Using this method, combined with an end to end optical 
model, we can determine the position of the BSM to within 
an equivalent of 0.5 arcsec on the sky, much better than 
the absolute pointing error of the Herschel satellite which 
is estimated as 2 arcsec (Pilbratt et al. [2010). 

4. Photometer bolometer time constants 

The basic time response of the photometer detector chan- 
nels is dominated by the low pass filtering in the signal 
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conditioning electronics which has been accurately mod- 
elled and verified during ground testing. However, the ther- 
mal response of the bolometers themselves has a second 
order effect that varies with bolometer temperature and it 
is important that this is understood as correcting for it is 
critical in both the astrometric and photometric accuracy 
achievable in SPIRE maps. The first assumption has been 
that the bolometer thermal response is represented by an 
'RC type response characterised by a single time constant. 
Investigation into the presence of any slower bolometer time 
constants on scan map and spectrometer time domain data 
has so far revealed no evidence for a slow response to op- 
tical power and only the correction for a first order fast 
time response is currently being performed by the pipeline 
processing. 

High energy particles hitting the bolometers cause sig- 
nal glitches as the energy is absorbed and the temperature 
of the bolometer increases. The response to these glitches 
gives a direct measurement of the bolometer and electron- 
ics time constants. To measure the photometer bolometer 
time-constant the data from a large number of photometer 
glitches have been collected and compared to the theoreti- 
cal impulse response modelled as a delta-function convolved 
with the photometer electronics and bolometer transfer 
functions. Applying this process to sets of individual (single 
event) glitches, reveals that the photometer glitch response 
for most photometer bolometers can be best described by a 
two-component transfer function, with a fast time constant 
T\ « 6 ms and a second slow component with a time con- 
stant of t 2 «105 ms and amplitude about 6.6 % (Fig. 1). A 
second type of common, concurrent glitches are also seen 
which are most likely due to interaction between particles 
and the bolometer array substrate. Applying the same pro- 
cess to these reveals that the common glitch response is best 
described by a first-order-only bolometer transfer function 
with a time constant of about 6 ms. This disparity in the 
system response between single-event and common glitches, 
and how this relates to the absence of any slow response to 
optical power, is not yet understood and is currently being 
investigated. The frequency response correction of both the 
electronics and the bolometers amounts to no more than 
5% of the signal for standard scan map data and given our 
knowledge of the electronics and bolometers any uncertain- 
ties in this correction can, again, be safely ignored in the 
overall calibration uncertainty. 



5. Photometer photometric calibration 

The conversion from voltage to flux density in the photome- 
ter is based on the assumption that the small signal re- 
sponsivity can be related to the bolometer operating point 
voltage nearly linearly i.e. that we can write the inverse 
responsivity in the form 



PSWB3 
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dV 
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V - K, 
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where S is the flux density, Ki,K 2 ,K 3 are constants and 
V is the measured voltage across the bolometer. K\ is in 
units of Jy/V, K 2 is in Jy and K 3 is in V. Integrating this 
equation gives 



S = K 1 (V-V )+K 2 \n 
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Fig. 1. The average response of a photometer bolometer to 
an ionising radiation particle derived from the co-addition 
of many glitches (crosses) compared to the modelled im- 
pulse response of the electronics and bolometer with (red) 
and without (blue) a slow time response component (see 
text). 



(2) 



where now Vq is another arbitrary constant that we chose to 
represent the bolometer voltage measured under a standard 
observing condition on a dark region of the sky. The vari- 
ous parameters have been determined from measurements 
in-flight using PCAL to measure the relative responsivity 
while viewing sources of different brightnesses in order to 
change the bolometer voltage. 

The calibration curves were measured for each detector 
in two steps. First, the change in voltage AVpcal observed 
for PCAL flashes were measured when the telescope was 
pointed at regions of the sky with different flux densities 
ranging from dark sky to the peak of Sgr A (for the nom- 
inal voltage settings) or Sgr B2 (for the lower gain bright 
source mode settings). The plot of 1/ AVpcal versus V 
gives the unsealed version of Eq. [TJ an example is shown 
in Fig. [2] where the telescope was pointed at the dark sky 
and different locations in the region of SgrA to vary the 
flux onto the detectors. The calibration curve itself is then 
scaled using data from observations in which Neptune (or 
a source calibrated using Neptune) is scanned across each 
detector in each array. The feedhorns and bolometers are 
not inherently sensitive to polarisation and any instrinsic 
polarisation in the calibration targets will have negligible 
effect on the overall photometric calibration. 

The initial absolute calibration scale was set using the 
asteroid Ceres and a model of its submillimctre flux den- 
sity based on the Standard Thermal Model (Lebofsky et 
al 11986ft . We have since observed Neptune and correction 
factors between the output of the photometer pipeline and 
the model of Neptune based on the work of Moreno (1998, 
I2TJTUJ) have been established (see Griffin et al. [2"0"T0]) . We 
discuss the accuracy of the Moreno model further in Sect. 
[7] and the accuracy of the calibration constants in Sect. [TU] 
of the present paper. 
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The fact that we do not require SCAL to null the tele- 
scope emission is a consequence of the lower total emission 
from the telescope and straylight compared to the expected 
values used in the initial design of the SPIRE instrument. 
Given this, and that using SCAL adds to the photon noise 
in the measurement, we have decided not to use the SCAL 
source in routine observations. An additional benefit of this 
mode of operation is that, with SCAL off it, and the rest of 
the instrument, are at a temperature between 4.5 and 5 K 
and the thermal emission from these components is limited 
to frequencies only detectable in the SLW band. 

The standard calibration for SPIRE spectra is based on 
a point source. An alternative calibration can be derived 
using the telescope which is appropriate for a source that 
entirely fills the detector field of view. A correction based on 
measurements of the spectrometer beam shape is required 
to convert from one to the other (see Sect. [7} . It should also 
be noted that instrument dependent spectral features in the 
passband of the spectrometer will change as a function of 
source extent. 



Fig. 2. An illustration of the basic method of photometric 
calibration. Here we plot the response to the internal cal- 
ibrator (1/AVpcal) versus the bolometer operating point 
voltage (V) as measured using flux densities ranging from 
dark sky to the peak of SgrA by pointing at different loca- 
tions in the region of SgrA. 



6. Spectrometer pipeline photometric calibration 

The calibration of the spectrometer follows a different 
method compared to the photometer. The signal that is 
measured by the spectrometer detectors is not a direct mea- 
surement of the flux density integrated over the passband 
as in the photometer but rather the Fourier component of 
the spectral content. Therefore Eq. [1] is not directly ap- 
plicable but an analogous equation can be used to correct 
for any non-linearity between absorbed power and bolome- 
ter voltage before transforming into the frequency domain. 
The parameters in this scheme are derived from a model re- 
sponse of the bolometer and have different values depending 
on the bolometer bias that is set. 

Once a linearised timeline in volts has been obtained the 
signal versus optical path difference is calculated using the 
mechanism position and further corrections for phase error 
are made (see Fulton et al. I2008[) . The source, telescope, 
SCAL and the instrument self emission are always mea- 
sured together and it is necessary to subtract a reference 
interferogram taken on dark sky to obtain the source on 
its own. This difference interferogram is transformed into 
spectral space and finally converted to flux density using 
a relative spectral response function (RSRF). The RSRF 
is derived by taking the interferogram of an astronomical 
source with a well modelled continuum, subtracting the ref- 
erence interferogram, transforming this into spectral space 
and dividing by a model of the source spectrum. In the 
present pipeline we use the asteroid Vesta as the calibra- 
tion source with a continuum model provided by T. Mueller 
(Mueller and Lagerros. l2002l> . 

Interferograms taken on Vesta, Neptune and Uranus 
showed that, with the SCAL source off, the signal at the 
central peak does not saturate, or at most only a few sam- 
ples are saturated, once the detector bias is correctly set. 



7. Testing the Astronomical Models 

To test the accuracy of the models used for both the spec- 
trometer and photometer we have derived the calibration 
using an alternative method where we use the telescope it- 
self plus knowledge of the instrument throughput (see Sect. 
IS]). As discussed in Sect. [2] the telescope is represented by a 
blackbody with a well known relative dependence between 
emissivity and wavelength even if the absolute overall emis- 
sivity may have an uncertainty up to 30% (see Fischer et al. 
2004). We can therefore generate an RSRF from the mea- 
sured telescope spectrum on dark sky and use this to cali- 
brate the spectra obtained on Vesta, Neptune and Uranus 
and compare to the models to check for self consistency. 
The results together with continuum models are shown in 
Fig. 131 The basic comparison is with the Neptune model 
of Moreno which shows an excellent agreement across all 
wavelengths giving confidence in both the relative and ab- 
solute level of the telescope emissivity. The Uranus model 
shown is also that of Moreno, but here we have increased 
the continuum by 10 Jy at all frequencies to give a better fit 
to the spectrum calibrated from the telescope. Reapplying 
a calibration file derived from this model of Uranus to 
Neptune shows that the Neptune measurement agrees with 
the model to within 5-10% across the full SPIRE band. The 
Vesta model as provided by Mueller was shown to be 15- 
20% lower than the measurement calibrated using Neptune 
and so we have increased the model spectrum that is used 
in the derivation of the spectrometer pipeline calibration by 
1.2. We discuss the accuracy of the pipeline derived spectra 
further in Sect. [TUJ 

8. Beam and astrometric calibration 

8.1. Photometer 

The measurement of the beam size of the photometer has 
been discussed in Sect. [2) The average beam derived for the 
three arrays has a near Gaussian core with a width of 18.1, 
25.2 and 36.6 arcsec respectively all with an uncertainty 
of ±5%. However, the full PSF outside of this core has 
much structure due to the diffraction from the secondary 
mirror and its support structure and it is recommended to 
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Fig. 3. The spectra of three solar system objects observed 
using the SPIRE Fourier transform spectrometer. From 
bottom to top we plot Vesta, Neptune and Uranus. The con- 
version to flux density was achieved using the Herschel tele- 
scope as a calibration standard as described in the text. We 
plot the model spectra of the three sources from respec- 
tively M oreno (J19981 I2010[) for Uranus and Neptune and 
Mueller (2002} for Vesta. The Uranus model spectrum has 
been increased by 10 Jy at all wavelengths and the Vesta 
model has been multiplied by a factor of 1.2. 



use equivalent beam areas of 501, 944 and 1924 arcsec 2 re- 
spectively to calculate surface brightness when observing 
extended sources (SPIRE Observer's Manual I20T0]) . These 
values are quoted as the mean for all detectors in an array, 
there is some minor variation across the arrays as demon- 
strated in Fig. [4] which shows the centroid positions and 
the relative beam sizes from a point source scanned across 
the PSW array. The variation in beam size has been found 
to vary by no more than 10% for the majority of detectors 
across any of the three arrays. In the example shown here 
the maximum and minimum beam sizes of the nominally 
operational detectors are 17.3 and 19.5 arcsec. The ap- 
parently large beamsizes in two of the detectors are caused 
by their very slow time response and these detectors are 
excluded from the map making. Measurement of the beam 
width versus wavelength within the bands is not possible 
in flight but characterisation was performed during ground 
test campaigns (Ferlet et al. l2008[) . 

The measured positions of the detectors in each ar- 
ray are folded into the map reconstruction algorithms 
together with the knowledge of the spacecraft pointing. 
Comparison between fields with sources seen at other wave- 
lengths (mostly from Spitzer and radio catalogues) shows 
that, with the present knowledge and algorithms, the typi- 
cal astrometric accuracy of SPIRE maps processed with the 
current version of the Herschel processing pipeline is around 
4 arcsec. The cause of part of the inaccuracy has already 
been identified and subsequent versions of the pipeline have 
an astrometric accuracy of 2 arcsec or better. This should 
be compared to the telescope relative pointing error of 0.2 
arcsec 1-sigma over 1 minute and the spatial relative point- 
ing error of 1.5 arcec (Pilbratt et al. |2010). 
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Fig. 4. A representation of the position (centres of the cir- 
cles) and beam size (radius of the circles) for the photome- 
ter PSW detector array. The missing circles indicate non- 
operational detectors. 
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Fig. 5. The measured spectrometer beam width as a func- 
tion of wavenumber for the central detectors in the two 
arrays. 



8.2. Spectrometer 

The shape and extent of the beam in the spectrometer has 
been subject to much study (Ferlet et al. 2008) both on 
the ground and in flight. The beam size versus frequency 
has been measured directly by taking medium resolution 
spectra on a point source of Neptune placed at different 
locations with respect to the central detectors by stepping 
the satellite position. The results are shown in Fig. \5\ The 
highly structured nature of the variation in beam size with 
frequency is expected from the multi-moded feedhorns used 
for the spectrometer arrays. The beams are only Gaussian 
at the low frequency end of each band and methods for 
conversion between calibration valid for point and extended 
sources require further investigation. 



9. Wavelength calibration 

9.1. Photometer 

The central wavelengths for the three photometer bands 
are derived from pre-flight measurements of the passbands 
of the detectors and the instrument optical filter chain 
(Spencer HEH SPIRE Observers Manual l2TJTtJ|) . As de- 
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scribed in Griffin et al. (2010) these are integrated over the 
waveband assuming a flat spectrum to give band centres at 
250, 352 and 504 /xm. The edges of the bands as defined at 
50% of the average in band transmission are 211-290 /xm 
(PSW), 297-405 urn (PMW), and 409-611 /im (PLW). The 
filter profiles themselves are available as standard calibra- 
tion products within the Herschel processing environment. 

9.2. Spectrometer 

The accuracy of the wavelength calibration was derived 
from line fits to 12 CO lines in the spectra of 5 sources: 
Orion Bar, CRL618, NGC7023, S106 and DR21. The space- 
craft velocity and the source velocity were subtracted from 
the fitted line positions to give the residual velocity off- 
set. Using this method we found a calibration accuracy 
within l/10th of a resolution element across the spectrom- 
eter band as defined by the FWHM of the Sine profile - i.e. 
1.207x0.04 cm -1 - and that the average line centre offset 
was 27 km s _1 . Much of this offset can be explained by the 
detailed physics of the interferometer (Spencer et al. [2010) 
which is now understood and we expect that all instrumen- 
tal offsets will be removed for all detectors across the field 
of view in future pipeline processing. The two spectrometer 
channels' transmission profiles fall rapidly at the edges and 
we define the trustworthy band edges for the central detec- 
tors of the two channels to be 14.9-33.0 cm -1 and 32.0-51.5 
cm . These edges vary somewhat depending on the posi- 
tion of detectors across the field of view but nowhere do 
they change by more than 0.2 cm . 

10. Discussion and conclusions 

The most significant part of calibrating an astronomical in- 
strument is the accurate conversion from measured signal to 
source flux density. In this paper we have given an overview 
of how this is currently carried out for the Herschel-SPIRE 
instrument. We have demonstrated an inter-comparison be- 
tween the models used for astronomical sources and the 
Herschel telescope emission spectrum that gives a final 
agreement between the Neptune model used as the pri- 
mary calibration standard and our measured spectrum of 
between 5-10%. 

The calibrator source model is only one component of 
the final uncertainty and we have tested the conversion pa- 
rameters for the photometer Volts to Jy using a Monte 
Carlo approach to determine the variance in the calibra- 
tion curve as a function of voltage. This method involved 
creating multiple calibration curves of 1/AVpcal versus V 
with addition random noise representative of the measure- 
ment errors and measuring the effect on the uncertainty 
in derived calibration coefficients. This translates into un- 
certainties that are < 0.1% for all detectors. This is much 
smaller than the other two sources of error, which are the 
uncertainty in the model for Neptune (5-10%) and the un- 
certainty in peak flux density measurements for Neptune 
derived from the maps after processing. Taking these fac- 
tors into account for data processed with the current pro- 
cessing software the overall photometric error can be as- 
sumed to be 15%. 

For the spectrometer photometric calibration a compar- 
ison of the pipeline derived Neptune spectrum based on 
Vesta to the Neptune model shows a typical variation of 
15-20% in the SSW band and 20-30% in the SLW band 



above 20 cm^ 1 and up to 50% below 20 cm -1 . It should be 
noted that variations in the temperature of the instrument 
itself between the measurement of the dark spectrum and 
the source can cause extra flux to be observed in the SLW 
band. The net result is an additive term in the spectrum 
below 25 cm -1 of the order of a few Jy which can leave 
significant positive or negative extraneous flux density lev- 
els in the source spectrum especially for faint objects. With 
careful interactive processing it is possible to remove this 
signature and obtain accurate source flux densities (Ivison 
et al. [2"0TUI Bocklee-Morvan et al. l2"UTo]) . Further investiga- 
tion on this aspect of the spectrometer calibration is ongo- 
ing and will be reported in future work. 

The calibration of the data from the SPIRE instrument 
processed through the current data reduction pipeline soft- 
ware is already of a standard more than sufficient to allow 
significant scientific conclusions to be drawn from the ob- 
servations. Future work will concentrate on improving the 
photometry using other astronomical sources and source 
models, including cross calibration with other facilities such 
as the Planck satellite. We will also work on improvements 
to the data processing tools, the astrometry of the maps 
through better attitude reconstruction and the wavelength 
calibration of the spectrometer through more sophisticated 
instrument performance models. 
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